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Abstract 

We demonstrate that the divergent terms in the OPE of a stress tensor and a surface operator of 
general shape cannot be constructed only from local geometric data depending only on the shape of 
the surface. We verify this holographically at d = 3 for Wilson line operators or equivalently the twist 
operator corresponding to computing the entanglement entropy using the Ryu-Takayanagi formula. We 
discuss possible implications of this result. 
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I. INTRODUCTION 


Global symmetries and the Ward identities that follow are extremely important in held theories. 
It gives stringent constraints to the allowed structure of correlation functions. In turn, given 
knowledge of how helds should transform under a symmetry transformation, this information has 
to be stored in the OPE coefficients between the global symmetry charge operators and the held 
operator under consideration. These OPE coefficients between charge operators and local operators 
have been studied extensively. To take the simplest example, in a 1+1 dimensional GET, the OPE 
between the stress tensor and a primary scalar operator takes the following form: 


T.z{z)0^{Q) 


A (Pa <9(Pa 

+-+ 


( 1 ) 


2;^ 2: 

where A is the conformal dimension of the operator (Pa- The hrst term stores data of how the 
scalar operator transforms under scaling transformation, while the second term describes how the 
scalar operator transforms under translation. One naturally wonders how the story generalizes to 
OPE’s of charges with surface operators. For starter, it is certainly the case that Ward identities 
still apply whenever there is a global symmetry. If a global symmetry acts locally, it is natural 
to expect that a local charge density operator can be dehned, and that the OPE coefficients - in 
this case coefficients of terms divergent in the inverse of the perpendicular distance of the charge 
density operator from the surface operator - should again carry information about how the surface 
operator transforms under the corresponding symmetry. These considerations should also have 
generalizations in studying higher symmetries [1] which necessarily involve symmetry transfor¬ 
mation of higher dimensional objects. OPE coefficients between the stress tensor and a surface 
operator is a very good starting point to study these questions, particularly so in a GET theory. 
The stress tensor exists in any GET and their expectation values are very strongly constrained by 
conformal symmetry, which allows us to make guesses of the form of these coefficients. Moreover, 
for co-dimension two operators wrapping spherical regions for example, the form of the expectation 
values of the stress tensor is constrained down to only a single coefficient, thus allowing us to make 
quantitative checks of our guesses. 

Go-dimension two operators are also very interesting for their connections to our understanding 
of entanglement entropy. The entanglement entropy and its cousins -the Renyi entropy- hold 
important information about the underlying held theory. It is known that the value of the entan¬ 
glement entropy can be understood as the expectation value of a twist operator that introduces 
branch cuts over the entangling surface surrounding the region under consideration [2-4]. An un¬ 
derstanding of these OPE coefficients of the stress tensor with the twist operators would allow us to 
compute at ease the entanglement entropy with entangling surfaces deformed arbitrarily from bet¬ 
ter understood ones, such as planar and spherical entangling surfaces. Ward identities of surface 
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operators, particularly those involving entanglement entropy, have been considered before[5, 6]. 
More recently there is a flurry of work towards understanding the entanglement entropy of more 
general shapes at least perturbatively from spherical/planar ones [7-12], Our work is the first 
attempt towards understanding these perturbative changes via OPE coefficients. 

Our strategy can be summarized as follows. We insert the stress tensor at a perpendicular 
distance e from the surface operator. As in [13, 14] we would like to focus on the ratio of the 
expectation value of the stress tensor to that of the surface operator alone. Simply from dimensional 
analysis, it is expected that the leading term of the OPE contains a 1/e'^ divergence as e —)■ 0, 
and that the form of the leading term is controlled by a single coefficient for a conformal field 
theory[15-17]. It is known that the leading term does not actually take part in controlling the 
transformation of the surface operator under scaling or translation[15]. This follows from Gauss 
law when we consider integrating the stress tensor over appropriate Gauss surface wrapping the 
surface operator. We are therefore most interested in the sub-leading terms in e. We make 
the working assumption that the the diverging sub-leading terms can in fact be constructed using 
local geometric data-such as extrinsic curvatures and torsion- of the surface on which the operator 
wraps. This approach was taken in [18] where the authors there attempted a systematic study 
of OPE’s of Wilson line operators in four dimensions, which is the classic example of extended 
operators. For simplicity, we focus on d = 3 theories, so that a co-dimension 2 operator is essentially 
a line operator, much like the Wilson line operator. We demonstrate however that our working 
assumption necessarily breaks down- that it is not possible to satisfy stress energy conservation if 
we insist upon constructing the sub-leading terms out of local geometric quantities. We check our 
claims by comparing our results with a holographic calculation of the stress tensor in the presence 
of a surface operator by inserting probes in the bulk. The holographic calculation produces an 
expectation of the stress tensor that is consistent with stress energy conservation. We demonstrate 
however that it is not expressible in terms of local geometric data, giving support to our field 
theoretic analysis. 

An outline of our discussion is as follows: we begin a review of the Ward identities of extended 
objects in section II, and then work out a systematic expansion of the expectation value of the stress 
tensor in the vicinity of a line operator in 2-1-1 d as a power series expansion in the infinitesimal 
perpendicular distance e, assuming that the terms are local geometric covariants. We obtain 
constraints of the coefficients of these terms up to the first subleading order in e using tracelessness 
and stress energy conservation. In section III, we compare these constraints with the known 
results of a line operator wrapping a circle, and which we find complete agreement, although we 
demonstrate why this program would begin to fail in the next sub-leading order in the e expansion. 
In section IV, we compute the stress tensor induced by the presence of a line operator slightly 
deformed by a transverse wave from a straight line holographically, and find evidence that support 
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our claim of a non-local stress tensor expansion. We conclude our paper with a few final comments 
in section V. Some details of the calculations are relegated to the appendices. 


II. WARD IDENTITIES AND OPES 


As shown in [5], under diffeomorphism, the entanglement entropy which is a special instance of 
a co-dimension 2 operator insertion, changes as follows: 

SEE(AJ*g)-SEE(A,g) = /(-^^) P) 

where / is a time-independent spatial diffeomorphsim continuously connected to the identity, and 
A denotes the region within the entangling surface on which the twist operator is inserted. Since 


SEE{f{A),g) — SEE{A,f *g), 


( 3 ) 


where f * g denotes the pull back of the metric, one can replace a smooth deformation of the 
entangling surface by diffeomorphism. 

The above is equivalent to 


\ y/Q / 




( 4 ) 


where S is the entangling surface, and F is some general function which depends on on coordinates 
tangential to S, and the unit normal vector, normal to the entangling surface and in the radial 
direction of the codimension 2 plane. This is the statement that the twist operator inserted over 
S breaks energy-momentum conservation, and momentum is only unconserved in the orthogonal 
direction. 

The l.h.s. is basically just the expectation value of the stress tensor in the presence of the twist 
operator. We would like to extract F, in which case we would be able to deal with any arbitrary 
deformation of the entangling surface, given the expectation value of the stress tensor. 


A. Constructing expectation value of the stress tensor relating to OPE 

From the discussion above, it is clear that at the end of the day, it is the UV divergent terms 
in the stress tensor expectation value that carries information about the functions F. Therefore, 
to extract this function for more general entangling surfaces, we would like to expand the stress 
tensor in powers of |e|. The expansion coefficients are then related to OPE coefficients of stress 
tensor with the surface operator inserted at the entangling surface. 
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To simplify the discussion, we will focus on d = 3, so that the entangling surface would be one 
dimensional. 

The expansion of the stress tensor should take the following form, 

T^y{x) = — H-^ . (5) 

We assume the OPE coefficients can be expressed in terms of local geometry data. Thus, for 
the leading term the most general form should be as follows: 

~ + 0-2tfj,tu + dsfluhty, ( 6 ) 

where t^{x) is the unit tangent vector on the entangling surface evaluated at the point x, and h is 
again the unit normal vector dehned by 

(x - x)^ 

rif, = - - 

e 

Some useful identities regarding the local geometrical quantities are relegated into the appendix 

A. 

Similarly, by using the dimensional analysis and parity symmetry, we can build out of the 
local geometric data as follows 

Al^ = K^ibiSf,^ + b2tf,U + b^h^n^) + bi{n^nlK^ + h^nlK^). ( 8 ) 

We note that the rule is that h and its derivatives thereof are allowed to appear an odd number 
of times at each order of the expansion, whereas and k must appear an even number of times, 
to preserve the symmetry of the problem. 

To continue, we also have 

^lu =codaK^n\^ty) + {ci6^^ + C2t^U + C2,n\nl) + 

+CiK\nlnlK^ + + c^n^nlK^K^ 

+ {K°'Y {cibfjiu + C2^ijty + c^hfj,hiy) . (9) 

We note that these parameters, {a^, bi, q, ■ ■ ■ } are not independent, since the stress tensor is 
conserved anywhere away from the entangling surface, and moreover in a conformal theory it has 
to be traceless. We can reduce the number of independent coefficients using these constraints, 
order by order in perturbation theory. 
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Plug in the expansion of one can obtain 


+ ■ ■ ■) + + \^\9^^liu + ■ ■ ■) • (10) 


Furthermore, 


and that 


AC 10 I “S"iP ( 1 -|<'|A'‘)) 

^ “ “U-|£|a-> + H ^ 


( 11 ) 


=(0cA'‘ - A'“nS0„n‘'’)(6i + " iC"') 


+a-‘('!#(2 


-) + fr 


nlK°‘ 


1 - leliPi' 1 - leliPi- 


+64 


^(2 - A II^J + ^« - - t.s-nl{K<^Y) 

\e\ i — |e| A ^ |e| 


( 12 ) 


where as dehned in (A7). 

Now collecting all the resnlts and reqnire energy-momentum conservation = 0, at the 

leading order of e-expansion of (10), i.e., (P(|e|“^), we have 


3oi -|- 203 — 0 . 


(13) 


To hrst snb-leading order, there are two independent terms whose coefficients we set to zero: 

— ti^K^ (3&1 bs dbi 03 ) -|- n^iF°^(a 2 + ^ 1 ) = 0 , (14) 

giving 

36i -|- &3 + 364 -|- 03 = 0, 02 + = 0. (15) 

Including constraints from tracelessness, we have in addition 

3ai -|- 02 T ®3 = 0, 3&1 -|- &2 T ^3 -|- 264 = 0. (16) 

In the next section we will check these relations by the known resnlt of a spherical twist operator, 
and then consider the more general snrface operators. 
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III. CHECKING THE CIRCLE IN A CONFORMAL FIELD THEORY 


In a conformal theory, we can determine the stress tensor in the presence of a spherical twist 
operator inserted along the spherical entangling surface. That is simply obtained by conformal 
transformation from the space H^-i X via 

—— =exp(-a), (17) 

ri UJ 

where u = r + it are the radial and Euclidean time coordinates respectively in the flat frame, 
and a = u + ir are the hyperbolic radial coordinate in Hd-i and Euclidean time coordinate of 
respectively. 

In the H^-i x frame, purely from symmetry one can deduce that the stress tensor is diagonal, 
taking the following form: 

TTT = —id—l)P, TAB = PgABi (18) 

where qab denotes the metric of H^-i 


Following through the conformal transformation, the stress tensor in the flat frame takes the 
following form: 


^ mYp 

{2RYP 

4|i?2_^2|d 


Ttr= i 


2R 
_ ^2 


2(2 -d)- 

(2(2-rf)+rf(p^ + 


dP fR^-dj^ R^-u^ 


R2 _ ^2 R2_ ^2 


/v» /v» /v» /v» 

^+4(5,,-^) 


These expressions can be expanded in e, where we take 

00 = R + {6r + idt), |e|^ = (5r^ + (22) 

from which we can extract the values of {ai,bi} and so on dehned in the previous section at least 
for d = 3. 




This allows us to evaluate the expected value of the stress tensor. In particular, 


+ = + (24) 

This can be compared with the corresponding expansion of the stress tensor 

(25) 

P / Sr‘^\ 

Trr= ^ - 26t^) - 3Sr—j +■■■ (26) 

Comparing the two expressions, we get 

ai = —2P, 03 = 3P = 02 (27) 

3 q P 

h = -3P, b, = -P = b,, &2 = ^. (28) 


We note that all the constraints (13) (15) (16) following from tracelessness and momentnm con¬ 
servation are satished as expected. We note also that these coefficients are over-determined and 
so provides a non-trivial check of the current construction. 


A. Hitting a rock— Failing Ward identity by local OPE for general entangling surfaces 

Despite that we have some success applying the ansatz to spherical entangling surface, challenge 
comes immediately when we try to deal with less special snrface S. The stress tensor from the 
ansatz is not conserved, nor does it provide the correct entanglement entropy pertnrbatively. 


1. Evidence 1 

In IIA we consider the energy-momentum conservation np to (!l(|e|“^ and it yields the con¬ 
straints (13)-(15). Now we will see that the conservation cannot hold at the order 0{\e\~‘^). At 
this order, the terms contribnted to are from the order (!l(|e|°) terms of 

/TAl, - rA'‘Al + \,\d“Al . (29) 

As we can see from (12), the term at this order from is , 

= (6i + b2)s,d^K^ + 0{K-f 
8 


( 30 ) 



where the hrst term is proportional to ti, and is of the hrst order of and the other terms 

involve but not d^K^. On the other hand, from the dehning equation (9) the terms in 

involving d^jK^ are instead in the form of which cannot compensate the 

term in (30). 

We note that in principle \e\d ^can also provide terms proportional to daK^t^. However 
as we can show explicitly that 

|e|9MjJa(|,|0) = 2c^K\K-nl - K^nl) + {K^)\c,nl + c^KK^) 

+ C4i^^ - nlK^ - Us'^nliK^f] + - {K^f]nl 

+ 2c'iK\K^nl - K^nl) - 2d[K\K^nl - K\l) 

+ {K^n4nl + c''n:Kn , ( 31 ) 

where we have used the identities in appendix A. We see no term in the form of which 

thus cannot be used to cancel the term in (30). Moreover, after detailed check we hnd that the 
other terms not involving will all be cancelled out exactly at this order. 

As a result, the conservation of stress tensor at this order at least requires 

6i + 62 = 0 , ( 32 ) 

which we know is inconsistent with 61,62 obtained from the case of spherical entangling surface 
(27). In other words, the stress tensor that follows from our ansatz is not conserved unless the 
entangling surface satishes d„K^ = 0 . 


2. Evidence 2 

In [11], it was shown that the entanglement entropy due to the deformation of a circular en¬ 
tangling surface (at d = 3) does not receive linear correction from the deformation. Instead, the 
quadratic correction is non-zero and takes the form 

ASee ^ - l){al + bl), (33) 

for the deformation in the form of 


Sr e{an cos n9 + bn cos n9) (34) 

where (r, 9) is the polar coordinate parametrizing the disk region bounded by the entangling circle. 
Similar contributions can be done for the deformation of a flat entangling surface. As we will 
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show from the holographic calculations in the next section, the result is 


SS 


EE 




+ bl), 


(35) 


for the deformation taking the form 

Sy = e{akCos kx + bksinkx) (36) 

to an entangling surface aligned along the x direction. 

We note that if this result is obtainable from the OPE of the stress tensor with the twist line 
operator, it has to come from as the linear order correction from vanishes. Namely, con¬ 
sider an entangling surface only slightly deformed from a flat entangling surface. The deformation 
is given precisely as a wave above. In this case, the extrinsic curvature of this surface must take 
the form 

K €k‘^{ak cos kx + bk sin kx). (37) 

Now if we further perturb by another Sy of exactly the same waveform, then to quadratic order in 
e, the only term that can contribute must be terms involving d^K ~ ek^{—aksinkx + bkcoskx). 
The hnal form of the surface integral of the stress tensor would have the term 

I ~ J dxd,K\^ ■ t) + 0{e ^), (38) 

which unfortunately vanishes at the order of 0{e^) because ■ t) ~ 0{e^). However, there is 
no other 0{e^) term from J , which is local and dependent only on the geometry of the 

entangling surface. 

This appears to be strong evidence that the OPE coefficients of the stress tensor must contain 
non-local terms, in addition to the terms we have found, despite the fact that the explicit results 
of a spherical entangling surface satishes all the constraints following from our local construction 
of the divergent terms. 


IV. HOLOGRAPHIC CALCULATION 

In this section we will re-consider the issues discussed in III in the holographic dual. Especially 
we will derive the holographic stress tensor in the presence of the entangling surface, i.e., the twist 
line operator at d=3, and find that it is conserved. The result implies that the OPE of the stress 
tensor with the twist line operator at d=3 cannot be in the local expression of geometric quantities 
as extrinsic curvature and its derivatives. 
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A. Planar entangling surface perturbed 


We first consider the variation of the holographic entanglement entropy due to the deformation 
of the flat entangling surface at d = 3. This is an analogy to what has been done in [11] for the 
deformation of the spherical entangling surface. 

We choose the following metric for the dual AdS 4 . 

ds^ = -17 {dt^ + dz^ + dx^ + dy^) (39) 

The entangling surface lies at y = 0 (also t = 0) and we consider the following perturbation 

5yQ = ee^^^ (40) 


for very small e. 

We can take the following ansatz for the corresponding minimal surface 

y{x,z) =yi{z)6yo (41) 

by also imposing the following boundary conditions to conform to (40) at z = 0, i.e., 

l/i(0) = 1, yi{oo) = 0 . (42) 


The minimal surface is obtained by minimizing the area 

v/(l + \d^y\^) (1 + \d^y\‘^) - \d^yd,y\^ 


A = 


dxdz, 


( 43 ) 


where we include the complex conjugate to make the area real. The change of the area of the 
minimal surface due to the deformation (42) is then given by 


AA = ^ [ dxdz 

/ 2z^ 


(44) 


from which we get the equation of motion for yi{z) 


zdlyi - 2a,j/i - h:‘zyi = 0 


( 45 ) 


Imposing the boundary condition (42), we obtain the special solution (assuming /c > 0 for 
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simplicity), 


yi(z) = (kz + l)e-'“. 


( 46 ) 


Substituting (46) into (44), we get the change of the holographic entanglement entropy 




EE 


AA = e" 


2/i 


j dx = e^ “ y + j dx . (47) 


This result is very similar to (44) in [11], Moreover, it is nonzero in contrast to the null result from 
the OPE of stress tensor with twist operator. 


B. Holographic computation of the stress tensor 

In the previous consideration we can obtain the change of the holographic entanglement entropy, 
bnt did not obtain the corresponding holographic stress tensor in the presence of the twist operator, 
which is the holographic dnal to the Ryu-Takayanagi (RT) action of the minimal surface. Now we 
would like to derive such a stress tensor holographically and check its conservation directly. The 
way to do it is to introduce a (probing) boundary source metric variation, which then backreacts to 
the bulk metric. We then evaluate the total on-shell bulk action, including the RT one [19] which 
represents the insertion of the entangling surface on the dual field theory side. That is, the total 
bulk action is Stotai ■= nSgr + (n — 1)See where Sgr is the bnlk gravity action and See the RT 
one. Taking the variation of the on-shell action w.r.t. the bonndary source metric variation and 
then setting the sonrce to zero, we will obtain the stress tensor in the presence of the entangling 
snrface inserted as a twist operator, i.e., the OPE of stress tensor with twist operator. 

We will not inclnde the backreaction of the surface operator on the background, presnmably 
because by taking the replica index n to one the dimension of the twist operator approaches zero, 
or eqnivalently in the holographic bulk its coupling to the backgronnd geometry is suppressed 
by n — 1. Since the extremal snrface satishes the equations of motion following from the Ryu- 
Takayanagi action, to linear order in the metric pertnrbation, we need only to evaluate the action 
on-shell neither backreacting the metric perturbation nor the extremal surface. 

We choose the Poincare coordinates for the AdS space, and the deformed metric due to the 
boundary sonrce metric perturbation takes the following form, 

^ [(1 + hzz)dz^ + {rjgy + hgy)dx^dx'' + 2hzgdzdxg\ (48) 
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We are following gauge choice in [20], in which 


where 


and 


^d-2 


hij Cd J d X Pijab dGab 


hzi Cd 


d 


^d-3 


d'^x 


hzz ^d 


d — 1 
d 


fd-i 

^d-2 


Biab dGab 


d-1 


f = , 


d X -y^Cab dGab , 


_ j-'tJ- 


Cd d f d\ ’ Biab ^diJj}^{di)Bj]^abi Gab JijiX) 

( 2 ) 

7 /-\ _ r „ 1 

d0 H.U 2“ 


1 Pijab i-j i^diadjl) “h dil) 6 ja) •> 


“ ^2 ’ ' 2 


(49) 

(50) 

(51) 

(52) 


and 5G^^ is the (probing) source metric variation at the boundary. We note that the expectation 
value of the boundary stress tensor is given by 


(T^")s = 


SSi 


total I 


ss 


EE I 


5G 


fiiy 


SG 


(53) 


flU 


In the above, ~ means that we take n —)■ 1 limit to suppress the backreaction to the bulk metric. 


We will focus on the case in which the entangling surface is almost planar, inserted at t := xi = 0 
and y \= Xd = ^- In this case the extremal surface is falling straight down, described again by 
7 / = 0 and extended along z,X 2 , - ■ ■ Xd-i- 


In the following we will consider only d = 3. We choose the static gauge such that the induced 
metric on extremal surface takes the form (up to the terms possibly contributing to 0(e)^ part of 
^See), 


kabda'^da’^ = ^ 


(1 + hzz + {dzyf + ^hzydzy) dz^ + 2h^zdxdz 


+ (1 + hxx + (Xy)"^ + ‘2hxydxy') dx"^ 


+ ■ • • , 


(54) 


where we have included derivatives of y w.r.t. the world-volume coordinates for latter convenience. 
For the strictly planar case, 7 / = 0 and therefore all terms involving y actually vanishes. 
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To the lowest order in e, i.e., 0(1), the stress tensor we get is that for a straight twist operator 




2 

StT (T2 + y'2)3/2 






(55) 


which by itself is conserved and traceless. 

The shape dependence of the stress tensor appears at the next leading order. Let us first 
consider the case in which the only nonvanishing component is SGxy The 0(e) order contribution 
to the (change of the) stress tensor reads {d = 3, x = x — x') 




5 


6G 


dxdz 


a /1 + hzz + 2iek[hxy + kz {h^y + + 0 ( e )2 


= e 


xy 

dxdz- 


SGxv~^^ 


ze 


ikx—kz 


TT^ + y'"^ + Z^y 


z + 1) + x"^ — 7/^ + z^^ 


-\-k ( t' + + y' + 


t'‘^{3kx + 2ikz + 2i) + 3kx'^ + 2ix^{kz + 1) 


-\-3kx [z^—y'^] + 2i{kz + 1) (+ 2 : 


/2 


(56) 


In the above we have set y{x,z) of the extremal surface into the form of (41) and (46). The 
integration over x in (56) can be done but we only manage to do the integration after an 
expansion in momentum k. After some tedious calculation, the result takes the following form. 



ee 


ikx' 


3it'‘^k 

StT (t'2 + |/'2)5/2 


i {3d^ + 2y'^) e 

IbTT 


Stt 


+ 0(k’‘) 


(57) 


The series in k above can also be understood as OPE of stress tensor with the twist operator 
since 5 = (f'^ + is essentially the minimal distance from the point where locates to the 

entangling surface {y = t = 0), i.e., the location of the twist line operator. 

The other components of stress tensor in the presence of the twist operator can be worked out 
similarly, and the detailed forms are put in the appendix B. Based on the detailed form of the 
stress tensor listed in appendix B order by order in /c-expansion, it is straightforward to show that 
the stress tensor is traceless and conserved, i.e., 

= 0 , = 0 . (58) 


Note that unlike the tracelessness, the conservation of is not easy to see without perform¬ 

ing the integration along the extremal surface. 

The particular form of the OPE coefficients seems to rule out the possibility that they be 
expressed covariantly in terms of local geometric quantities. For example, we can consider the 
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term in (Tyy^^^^) (see (B3)) and denote the coefficient as . Note that the extrinsic 
curvature K in this case is and none of ty have any components proportional to 

k'^. In other words, a^y should contain one K^. On the other hand, is also proportional to y^, 
or rather (n^)^ and does not appear to follow from any 2-tensor built from and 


With one can recover the holographic entanglement entropy computed in sec IV A. 

More explicitly, we consider perturbation of the entangling surface given by 6y = (i.e., 

e —)■ e -|- e). With the following coordinate transformation 

y ^ y — z ^ z, x ^ x , (59) 


the entangling surface can be made as a plane {y = 0) again. The resulting hrst order perturbation 
to Gyu is 

5G^y = -ikeG^\ (60) 

with all other components of 5Gij being zero. This perturbation can also be understood as being 
generated by the vector ^ whose only nonvanishing component is ^y = —Note that with 
an extra gauge transformation = ikye^^^, we can get to the metric in the normal coordinates 
used in [7] (which has Gxy = 0, Gyy oc yK^^). This vector ^ vanishes on the entangling surface 
{y = 0) and therefore the diffeomorphism it generates does not change the entanglement entropy. 
Substituting in (54) and (60), we get the leading order correction to entanglement entropy. 


ASee= / d^x'{Txy‘^^^\x'))6G:yix') 


= ee y d^x' 
1 


dxdz 


Qikxz 


2{kz + 1) (t'2 + - 7y'^ + z^) 


^ 2/4 


e {t'‘^ + x^-y'^ + z^) + 


^2f5 

Ak'^z{kz + 1) 


TT^/^ 


Akx—kz 


+ C>(£=) . 


(61) 


Note that y{x, z) in hzi only gives contribution of higher orders in e and therefore can be set to be 
zero. On the other hand, unlike hzi, hzz is not multiplied by e and therefore we need to expand 
y ^ 0{e) in / and Gxy The integration of y' and t' can be performed straightforwardly and the 
result is 


AS'rr = ee / dx' / dxdzk'^z 


—2ikx 


2{kz + l) 

TT (£2 -7 z'^)'^ 71 (£2 -|- z'^)'^ 


+ 


Akx—kz 


(62) 


We notice that hzz actually has no contribution to ^Aeb- The integration of x is essentially the 
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Fourier transform and we get 


ASee = ee/c 


= eek 


^OO / 

^ J (k‘^ + 

(1 + kjji)e~‘^^^ 


{kz + 1)^ 


dx' 


/r 


J dx' = 2ee ~ "Y J > 


(63) 


where /i is the IR cutoff in the bulk. The result exactly reproduces (47). To be more precise, we 
perform the replacement e —)■ e + e in (47) and hnd the terms of the order of 0(ei) agree. 


V. CONCLUSION 

In dimensions higher than two, a co-dimension two twist operator is a surface operator and 
the precise formulation of its OPE with some local conserved charge density such as stress tensor 
have not been fully explored. As a starter, in this paper we study the OPE of stress tensor with 
a twist line operator in 3 dimensions. We try an ansatz for such OPE coefficients constructed out 
of local geometric data encoding the shape of the line operator. It turns out that the resulting 
stress tensor in the presence of the twist operator can not be conserved for general (wavy) shapes 
of line operator except for highly symmetric one such as planar or spherical line operators. This is 
in contrast to the conserved holographic stress tensor, which we have constructed in this paper in 
the presence of the twist line operator based on RT formulation. Furthermore, we demonstrated 
that the stress tensor obtained holographically is not expressible as local geometric expansion, 
implying that the OPE of stress tensor with line twist operator cannot be expressed in terms 
of local geometric data only. It is possible that the Ward identities could be anomalous in the 
presence of surface operators, however, the holographic results suggest otherwise. 

The perturbative computation we perform in sec IV is done holographically. It would be inter¬ 
esting to see whether the same stress tensor and entanglement entropy can be obtained from the 
integrated correlation functions of the CFT side using the second order perturbation (two insertions 
of stress tensors, see e.g. [7]). However, it should be noted that there might be some subtleties in 
pushing the perturbative computation to second order. For example, it remains unclear [10] how 
to reproduce the shape dependence of the entanglement entropy (in four dimensions) using this 
approach. 

We should emphasize the failure of OPE in terms of local geometric data is peculiar in odd 
dimensions. In contrast, the universal part of the entanglement entropy in four dimensions is known 
to have local expressions [21, 22] in terms of (Graham-Witten) conformal anomalies encoding shape 
information of the twist operator. Thus, on hindsight, our results are not unexpected as there is no 
conformal anomaly in odd dimensions so that the universal part of the entanglement entropy has 
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no even dimensional analogue. Thus, it deserves further study of the issues raised in the current 
paper to get a clearer understanding of the results in both even and odd dimensions. 
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Appendix A: Useful Identities regarding the local geometrical quantities 

To obtain further expansions, we would need a few more building blocks. 

To begin with, we need a complete basis of vectors in the vicinity of x, let us also dehne the 
vector 

(-^ 1 ) 


One can readily check that = 1 and by construction it is normal to both and t 
We can evaluate the extrinsic curvatures, dehned as 


t^- 




(A2) 


where n“, a = {1, 2} denotes the two possible unit normals = n and = k respectively. 
One very useful identity is 


da 


t„ 


^X^^ 1 — n°‘data ’ 

where a is a convenient parametrization of the entangling curve, such that it’s trajectory in space- 
time can be denoted by {y^{a)}, and that 


= w 


da ’ 


(A4) 


where the normalization factor A/” = 1 if a is chosen to be some affine parameter. 
The extrinsic curvature is therefore given by 


= <(9X, 


(AS) 
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for affine parameter a. 

This identity is derived by requiring that as the insertion point is varied, x{ao) should vary 
accordingly so that 

= 0 (A6) 


remains true. 

There are a few useful identities: 


d^nl= 

d^K^= 

d,t^= 

dt^_ 

da~ 

d^nl= 


H A" - - ’V".). = = rvyjTfi ■ 

- nlnlt/s) + 

f)2 n 

77 ^ 3 = -S — 0 

'V’ ^ 

+ KV 

d^v’' 

^ s =n‘^K°-H 
Q^2 

= <ir“, 


(9(t2 

= nlOy 


-t ^_|^|^i- 


(A7) 

(A8) 

(A9) 

(AlO) 

(All) 

(A12) 

(A13) 


Appendix B: Holographic stress tensor in the presence of twist operator 


Follow the same procedure in deriving , one can derive the other components of the 

stress tensor in the similar way. The hnal results are listed below by /c-expansion and integrating 
out the worldvolume f dx ■ ■ ■, 


0(e)^ _ ^^ikx 


ee 




i (3(2 + 2!/2) e 


{r,yi''> = ee 


ikx 


StT IStj- (^2 _|_ StT 

3t {f - 47/2) 3t {f + 27/2) A;2 


+ 


tk^ 


StT + 7/2)A2 IGtT + 7/2)^A 47r(t2 + 7/2) 

t (9t2 + 87/2) ^4 


+ C>(A;'' 


(T„OI''> = ee 


Sy^k"^ 


+ 


yk^ 


6477 + 7/2)^/2 

37/ (7/2 - 4t2 

StT + 7/2)A2 ^^2 _|_ 4:7r{t‘^ + 7/2) 

7/(12t2 + 117/2) 


6477 + 7/2)^/2 


+ c>(r 


(Bl) 


(B2) 


(B3) 
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and 






ee 


ikx 




Syk"^ 


+ 


5yk* 


ee 


ikx 


Stt (t^ + IGtt (t^ + 647rA/t^ + y'^ 

3 (3t^|/ — 2y^) 3y (t^ + 2y‘^) k"^ yk^ 

Stt (t^ + y‘^y^‘^ IGtt (t^ + ^tt (t^ + |/^) 

y iyt^ + 6|/^) 


+ C>(A;® 


647r (t^ + y’^y^'^ 

3ityk 


ee 


ikx 


+ 0{k^ 


+ 


ityk^ 


Stt (t^ + y‘^y^‘^ IGtt (t^ + 


+ C>(A;® 


where we have dropped the primes of t' and y' for simplicity. 


(B4) 


(B5) 

(B6) 
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